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ABSTRACT 

We present HST WFPC2 I band imaging for a sample of 9 Hyperluminous Infrared 
Galaxies spanning a rcdshift range 0.45 < z < 1.34. Three of the sample have mor- 
phologies showing evidence for interactions, six are QSOs. Host galaxies in the QSOs 
are reliably detected out to z ~ 0.8. The detected QSO host galaxies have an el- 
liptical morphology with scalelengths spanning 6.5 < r e (Kpc) < 88 and absolute k 
corrected magnitudes spanning —24.5 < Mj < —25.2. There is no clear correlation 
between the IR power source and the optical morphology. None of the sources in the 
sample, including F15307+3252, show any evidence for gravitational lcnsing. We infer 
that the IR luminosities are thus real. Based on these results, and previous studies 
of HLIRGs, we conclude that this class of object is broadly consistent with being a 
simple extrapolation of the ULIRG population to higher luminosities; ULIRGs be- 
ing mainly violently interacting systems powered by starbursts and/or AGN. Only a 
small number of sources whose infrared luminosities exceed 10 13 Lq are intrinsically 
less luminous objects which have been boosted by gravitational lensing. 

Key words: infrared: galaxies - galaxies: active - galaxies: Seyfert - galaxies: star- 
burst - Quasars: general - gravitational lensing 



1 INTRODUCTION 



Serjeant ct al 1995 ; Eisenhardt et al 1996j ; |Green fc Rowan- 



Robinson 1996). These objects appeared to presage an en- 



One of the most important results from the Infrared Astro- 
nomical Satellite (IRAS) all sky surveys was the detection 
of a new class of galaxy where t he bulk of the bolometric 
emission lies in the infrared range ( Soifer et al 1994 ; Banders 
& Mirabel 1996). This population becomes the dominant ex- 
tragalactic population at luminosities above Lm > 10 11 Lq. 
At the extreme upper end of the IRAS galaxy population 
lie the Hyperluminous Infrar ed Galaxies (HLIRGs), those 
with Lir > 10 13 /i e r 5 2 L Q ( |Rowan- Robinson 2000] ). The 
first HLIRG to b e foun d, P09104+4109, was identified by 
Kleinmann et al fll988] ), with a far infrared luminosity of 



1.5 x 10 13 ft 50 2 Z/Q. In 1991, Rowan-Robinson et al identified 
F10214+4724 with z = 2.286 and a far infrared luminosity of 
3 x 10 h^£_L0. Later observations of this object revealed a 



huge m ass_ ofjgokcuh^j^ fc vandcri 



Bout 1991| ; |Solomon, Downes fc Radford 199$ ), a Seyfert 
emission sp ectrum (lElston et al 1§H ) , very high optical po- 
larisation ( Lawrence et al 1993 ), and strong evidence for 
lensing with a magnification o f about 10 at infrared wave- 
lengths (Graham fc Liu 1995; Broadhurst & Lehar 1995 



tirely new class of infrared galaxy. 

There are currently about 50 sources positively identi- 
fied as HLIRGs, of which 13 were identified in spectroscopic 
follow-up observatio ns of sources discovered in far infrared 
or sub- mm surveys (Rowan- Robinson 2000). Two of those 
13 IR-selected HLIRGs have been found to be lensed. By 
projecting population de nsities from the IRAS Faint Source 
Survey Rowan- Robinson ( 2000 ) estimates that there are be- 
tween 100 and 200 HLIRGs with 60/im flux > 200 mjy over 
the whole sky. 

The source and trigger of the IR emission in HLIRGs is 
currently the subject of considerable debate. ULIRGs ap- 
pear to be powered by s tarburst and/or quasar activity 
trigg ered by interactions (Leech et al 1994; |Clements et a] 
i99ej ), the HLIRGs may simply be the high lumi nosity tail of 



the ULIRG population. Rowan- Robinson (2000) argues that 



the majority of the emission at rest-wavelengths >50/im in 
HLIRGs is due to starburst activity, implying star formation 
rates > 1000MQyr _1 . If the rest-frame far infrared and sub- 
mm emission from HLIRGs is due to star formation, then 
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Table 1. Hypcrluminous Galaxies Observed By HST 



Name 


RA (2000) 


Dec 


z 


f 1 

•<814 




M 2 
I 


L 3 

IR, 


IR SED 7 


Opt Spec. 8 


Morphology 


F00235+1024 


00 26 06.7 


10 41 27.6 


0.58 


0.08 


18.73 


-23.6 


13.1 4 


63% Sb, 37% AGN 4 


NL 


Interacting 


F10026+4949 


10 05 52.5 


49 34 47.8 


1.12 


0.04 


19.61 


-24.3 


13.8 5 


AGN dominated 5 


Syl 


Interacting 


PG1148+549 


11 51 20.4 


54 37 32.8 


0.97 


2.07 


15.21 


-28.3 


13.7 


30% Sb, 70% AGN 


QSO 


Undisturbed 


LBQS1220+0939 


12 23 17.9 


09 23 07.3 


0.68 


0.43 


16.91 


-25.8 


13.5 5 




QSO 


Undisturbed 


F12358+1807 


12 38 20.2 


17 50 38.8 


0.45 


1.09 


15.90 


-25.9 


12.7 




BAL QSO 


Undisturbed 


F12509+3122 


12 53 17.6 


31 05 50.5 


0.78 


0.91 


16.10 


-26.9 


13.5 5 


AGN dominated 5 


QSO 


Undisturbed 


F14218+3845 


14 23 55.5 


38 31 51.3 


1.21 


0.07 


18.96 


-25.1 


13.1 4 


74% Sb, 26% AGN 4 


QSO 


Undisturbed 


F15307+3252 


15 32 44.1 


32 42 46.2 


0.93 


0.14 


18.20 


-25.2 


13.5 4 


32% Sb, 68% AGN 4 


Sy2 


Interacting 


PG1634+706 


16 34 28.8 


70 31 32.8 


1.33 


6.30 


14.00 


-30.4 


14.0 


17% Sb, 83% AGN 


QSO 


Undisturbed 



Coordinates and absolute F814 W band magnitudes/fluxes are taken from the HST images and are in the VEGAMAG HST flight filter 
system. 1 1 band flux in units of mjy. Subscript denotes wavelength in microns. 2 No applied k correction. 3 Logarithm of the 1 — 1000/xm 
IR luminosity, where necessary recomputed for Ho = 65 km s~ 1 Mpc - 1 and Qo = 1.0, and in units of bolometric Solar luminosities. 
4 Verma et al 2001. 5 Rowan- Robinson 2000. 6 Haas et al 1998. 7 IR Spectral Energy Distribution, computed between 1 — 1000/^m. 8 Optical 
spectral classification, taken from Rowan-Robinson 2000. 



the implied star formation rates would be the highest for 
any objects in the Universe. This would strongly suggest 
these galaxies are going through their maximal star forma- 
tion periods, implying that they are very young galaxies. 
A third possibility is that if the IR emission arises via some 
other mechanism (e.g. a transient IR luminous phase in QSO 
evolution) then these galaxies may represent an entirely dif- 
ferent class of object. 

In this paper we examine the morphologies of HLIRGs, 
their host galaxies and immediate environments in an at- 
tempt to discern the nature and origins of their IR emission. 
We present a sample of Hyperluminous Infrared Galaxies im- 
aged with the Wide Field Planetary Camera 2 on board the 
Hubble Space Telescope. Sample selection and observations, 
and data analysis are described in §§2 and 3. Results are 
presented in §4, including sample morphology, lensing prop- 
erties and descriptions of individual sources. Discussion of 
these results can be found in §5 and conclusions can be found 
in §6. 

Unless otherwise stated we have taken Ho = 65 km s -1 
Mpc -1 and fl = 1.0. 



terparts. LBQS1220+0939 and F12358+1807 are associated 
with known quasars with tub < 18.5, i.e. 10 times less bright 
in the rest frame UV than the PG quasars. F12358+1807 
lies slightly below the 10 13 Lq limit but is a member of the 
Broad Absorption Line quasar class. BAL QSOs are believed 
to exhibit excess FIR emission and F12358+1807 is the most 
l uminous object in its class. F15307+3252 is an IRAS source 
( |Cutri et al 1994 ) that has an IR luminosity in excess of 



1O"L . 

The HST data were taken in cycle 6 between March and 
July 1997 using the Wide Field Planetary Camera 2. The 
coordinates of each source were centred on the Planetary 
Camera CCD, selected for its superior pixel scale and full 
sampling of the PSF over the Wide Field Camera CCDs 
(0.046" pixel -1 as opposed to ~ 0.1" pixel -1 ). Exposures 
were taken using the F814W filter, which corresponds closely 
to the / band filter in the Cousins system. Each observation 
consisted of one short exposure (100 - 600s) to image the 
brightest parts of each object with minimal saturation of the 
detector, and two or three longer exposures (700 - 1600s), 
taken both to image the fainter regions of each source and 
to facilitate the subtraction of cosmic rays. 



2 OBSERVATIONS 
2.1 Sample 

The 9 objects in this sample all have a far infrared luminos- 
ity close to or exceeding 10 13 Lq that cannot be explained 
by non-thermal emission. The redshift range spans z = 0.44 
to z = 1.34, with a mean redshift of z = 0.9. Of the ob- 
jects in this sample, 7 were identified purely by IR selection 
techniques. 

Three of the sources, F10026+4949, F12509+3122 and 
F14218+3845 were discovered as part of a redshift survey 
of 3703 IR AS FSS objects to a flux limit of 200mJy ( p liver 
et al 1996]). PG1148+549 and PG1634+706 are IRAS de- 



tected PG QSOs with SEDs indicating an IR luminosity 
of Lir > 10 13 L Q . F00235+1024, LBQS1220+0939 and 
F12358+1807 were discovered as the result of a systematic 
program of statistical optical identification being carried out 
with the APM machine (McMahon et al 2001, in prepa- 
ration). This project uses robust statistical estimators to 
identify IRAS sources taking into account the IRAS error 
ellipse and optical magnitudes of all potential optical coun- 



3 DATA ANALYSIS 
3.1 Data Reduction 

All datasets were calibrated on arrival with the best avail- 
able reference files, using the IRAF task CALWP2, and com- 
bined into a single image using the IRAF task GCOMBINE. 
Statistical weights were assigned to different exposures 
based on their exposure times. The GCOMBINE rejection al- 
gorithm CCDCRREJ, which incorporates the Planetary Cam- 
era CCD characteristics, was used to remove cosmic rays by 
comparison between different exposures. Pixels flagged as 
saturated in the deeper exposures were replaced by their 
scaled unsaturated counterparts from one of the shallower 
exposures, with suitably modified statistical weights. Fol- 
lowing these steps the sky background was subtracted and 
warm pixels were removed by linear interpolation from their 
immediate neighbours. 

Photometric solutions were calculated using the Syn- 
phot package from the STScI, whic h is c onsistent with the 
solutions given by Holtzman et al (1995), and are given in 
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the VEGAMAG system. For the QSOs in the sample that 
were saturated in the central regions PSF fitting photom- 
etry was performed to accurately measure the true source 
magnitudes. Corrections were made for detector gain (for 
these observations the gain was 7e~ / DN) and aperture size. 
Absolute magnitudes (without a k correction term) were de- 
rived using the expression: 

5log[(2c/H )(l + z)(l - (1 + z)- 1/2 )] (1) 



M 



25 



We estimate that our derived relative magnitudes have an 
associated photometric error of approximately 5%. Prescrip- 
tions for converting magnitudes in the HST filter system 
(the 'flight' system) into more conventional filter systems 
depend in some way on the object colours. The F814W filter 
is a good match to the Cousins / band filter, the conversion 
factor is very small (< 0.05 magnitudes) and depends only 
slightly on the object colour. As the colours of our sample 
are not known we have not applied a conversion factor and 
our magnitudes, although referred to as / band magnitudes, 
remain in the WFPC2 flight filter system. The difference be- 
tween these magnitudes and the equivalent Cousins / band 
magnitudes is small (< 0.1 magnitudes) and thus compar- 
isons can readily be made between the two systems. Fluxes 
in Jy were obtained by multiplying the flux in ergs cm -2 
s _1 A -1 by the central wavelength of the F814W filter in 
angstroms, and then dividing by the central frequency in Hz. 
As such they do not include any bandpass colour correction. 

k corrections for the F814W filter were computed using 
the spectral synth esis package PEGASE and the as sociated 
galaxy temp lates (Fioc fc Rocca-Volmerange 1997 ). Previ- 
ous studies ( purace et al 1998| ) have shown that the optical 
emission from very luminous IRAS galaxies is dominated by 
emission from old stellar populations, rather than an unob- 
scured starburst. We therefore used models for evolved stel- 
lar populations to calculate k corrections. The galaxy tem- 
plate for a stellar population of age 10 9 years was generated 
using a Rana & Basu (1992) IMF and a star formation law 
r(t) = vg{t) where g(t) is the gas fraction. Published val- 
ues of k corrections can vary markedly, especially for bluer 
filters (e.g. Poggianti 1997). We estimate that the error in 
our computed values due to uncertainties in Initial Mass 
Functions, star formation laws, the age of the system and 
the assumed source spectrum is approximately 0.2 magni- 
tudes at z — 0.45, rising to 0.45 magnitudes at z — 1.33. 
No corrections have been made for reddening due to dust 
extinction. 



3.2 The Point Spread Function 

Host galaxies can only be resolved in QSOs if light from 
the central bright source, and from the host galaxy, can be 
separated. The most convenient way to accomplish this is 
by fitting a point source template to the QSO; the template 
being established from a set of observations of stellar Point 
Spread Functions (PSFs) that are a good match in filter, 
colour and detector position to the targets. For the objects in 
this sample no additional observations were made of nearby 
stars, and no suitable observed PSFs were available in the 
HST PSF Archive. Synthetic PSFs were there fore generated 
using the Tiny TIM v5.0 software ( Krist 1995 ), recently up- 
dated to include field dependent effects such as astigmatism, 
coma, and focus. 



Ten times oversampled PSFs were generated. These 
PSFs were then shifted by non-integer pixel distances, re- 
binned to normal resolution and convolved with the PC 
pixel scattering function. The agreement between these syn- 
thetic PSFs and observed (stellar) PSFs is generally excel- 
lent within radii of about 2". Beyond this radius scattered 
light in the WFPC2 optics can have a significant effect, an 
effect that is not modelled by the Tiny TIM software. The 
effect of any systematic uncertainties in the PSF have also 
not been examined, such effects are however likely to be 
extremely small. 

The resulting PSFs were then centroided and normal- 
ized to the image via a reduced \ 2 fit- It was found that all 
the images were saturated in the central regions to varying 
extents, mostly confined to the central 3x3 pixels and never 
beyond the central 5x5 pixels. In all cases the centralmost 
possible pixels were used, avoiding those that were saturated 
or that lay in the diffraction spikes. This method allows the 
image and PSF to be registered to within ~0.2 pixels (esti- 
mated by visual inspection), and also gives a starting value 
for the PSF normalisation. 



3.3 Profile Fitting 

In order to determine the morphology of the host galaxies 
in the sample, both de Vaucoleurs and exponential disk pro- 
files were fitted to the surface brightness profiles of the host 
galaxies. The de Vaucoleurs profile is of the form: 



J(r) 



j -7.67[(r/r e ) 1 /4„ 1] 



(2) 



and the exponential disk profile is of the form: 

I(r) = I e' {r/h) (3) 

One dimensional surface brightness profiles of the 
sources were extracted by fitting annuli spaced at 1 pixel in- 
tervals, fixing the ellipticity (at e — 0.1), semimajor axis and 
position angle. The PSF profile together with either a radial 
de Vaucoleurs or disk profile were then fitted to the source 
surface brightness profile, treating the PSF normalisation as 
a third free parameter. This helps to minimise as far as pos- 
sible the systematic errors introduced by subtracting off a 
'best guess' PSF and fitting a galaxy profile to the remaining 
light distribution. Fitting was carried out using Levenberg- 
Marquardt least-squares minimization in order to robustly 
derive the best fit parameters. Two dimensional isophote 
fitting, using the IRAF task ELLIPSE, was also used to in- 
vestigate the ellipticities and position angles for each host by 
iteratively fitting isophotes of constant surface brightness. 

In order to quantify whether a host galaxy was detected 
or not we imposed the condition that a PSF + galaxy to 
QSO fit must be better than a pure PSF fit to at least 
95% confidence, computed via an Ftest. We excluded pixels 
within the central 0.1" as pixels within this radius tend to 
be highly undersampled and in some cases saturated. We 
also excluded all pixels beyond 2" because of light scattered 
within the WFPC2 optics. 



4 RESULTS 

Coordinates, redshifts and measured F814W magnitudes for 
the sources can be found in Table [l], together with the opti- 
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Figure 1. HST F814W images of the 9 HLIRGs, showing the details of each source and scaled to show any interesting features. 



cal morphologies, IR luminosities and any previously cited 
best fit IR SED models. Table |^ summarises the host galaxy 
properties, and includes relative/absolute magnitudes, k cor- 
rections and host galaxy properties derived from profile fit- 
ting. In the sections where host galaxies are discussed, this 
refers to those objects listed in Table 2, thus including three 
sources which do not possess a central bright source. Discus- 
sion of PSF subtraction and/or host galaxy detection limits 
are not applicable to these objects. 



4.1 Morphologies and Lensing 



The HL 'T images showing the sources and immediate envi- 



source, F10026+4949 is small, shows no bright central 
source and slightly weaker signs of interaction. One object, 
F14218+3845, is a QSO with a relatively underluminous 
point source. The five remaining sources are all optically lu- 
minous QSOs. Although some of the QSOs show some very 
faint inhomogeneities, none of them show any definite signs 
of recent or ongoing interaction. 

None of the objects in the sample show any evidence 
for gravitational lensing. We investigated the possibility of 
lensing for all sources using the HST gravitational lens sim- 
ulator constructed by Ka van Ratnatunga as p a rt of the HST 
Medium Deep Survey ([Griffiths et al 1994 ; Ratnatung; 



Griffiths & Ostrander 1999). For non-trivial magnification 



ronments can be found in Figure |2j. Detailed images of each 
object can be found in Figure ^. 

The HST images show that the sample contains an 
interesting mixture of morphologies. Two of the sources, 
F00235+1024 and F15307+3252 show clear signs of strong 
interaction, features include obvious morphological dis- 
turbance, tails, and compact bright 'knots'. One further 



arcs of length comparable with the Einstein radius should 
be observed, so to first order we can set limits on the Ein- 
stein radius of 0.3" (~6 PC pixels). If any lensing is there- 
fore present it would have to be by low mass galaxies and so 
would make no appreciable difference to the derived optical 
and IR fluxes, unless the background source happened to be 
aligned with a caustic. The lack of multiple images around 
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Figure 2. HST F814W band images of the 9 Hyperluminous galaxies, showing the sources (circled) and immediate environments. 



any of the QSOs limits any possible lensing magnification to 
less than a factor of 2 in these objects, assuming an isother- 
mal sphere model. If no other signs of lensing are apparent 
at HST resolution however, such as tangential stretching, 
then it is extremely unlikely that lensing significantly con- 
tributes to the luminosity of any of the objects in the sample. 
Those sources whose morphologies are, on first impressions, 
plausibly consistent with some form of lensing scenario are 
discussed further in §4.5. 

For two of the three sources not dominated by a bright 
point source, F10026+4949 and F15307+3252, we were able 
to fit theoretical galaxy surface brightness profiles. For both 
these sources the galaxies were found to be extremely lu- 
minous ellipticals, the properties of which can be found in 
Table §. The morphology of F00235+1024 is so disturbed 
that no type of galaxy profile can be fitted to the source 
surface brightness profile. 



For the two PG quasars Schmidt & Green ( 1983 ) have 
measured fluxes at 440nm. Their results, compared to the 
measured F814W fluxes given in Table hi, suggest that the 



optical continuum for both these sources is extremely flat, 
with a ~ 0. 



4.2 Infrared Properties 

Table |l] lists the IR properties of the sample. Three of the 
sources (F00235+1024, F14218+3845 & F15307+3252) have 
been previously observed by the Infrared Space Observa- 
tory (ISO) as part of a larger sample, and their Spectral 
Energy Distributi ons modelled usin g a suite of AGN and 
starburst models (Verma et al 2001). The IR emission from 
PG16 34+706 has been previously modelled by Haas et al. 
using a series of greybody functions. The IR emission 



dl998Q 



from a further three sources (F10026+4949, F12509+3122 
& PG1148 +549) has been previously examined by Rowan- 
Robinson ( 2000[ ). For the final two sources in the sample 
(LBQS1220+0939 & F12358+1807) there exists insufficient 
photometry to model the power source behind the IR SED. 
Some constraints on the IR emission from F10026+4949 & 
F12509+3122 can also be drawn using the IRAS colours. A 
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'warm' infrared colour (i.e when the IRAS 25^im to IRAS Table 2. Host Galaxy Properties 

60/im flux ratio has a value of /2s //60 > 0-2) generall y 



indicates the presence of an AGN (de Grijp et al 1985). 
Both F10026+4949 (/25//60 = 0.665) and F12509+3122 
(/25//60 = 0.472) possess 'warm' colours, and have been 
previo usly cited as containing an AGN by Rowan- Robinson 
(pOOOh. For F12358+1807 and LBQS1220+0939 the IRAS 



fluxes and upper limits do not allow us to distinguish be- 
tween 'warm' and 'cool' IR colours, hence the power source 
is unconstrained. 

For two of the sources in the sample (PG1148+549 & 
PG1634+706) we present new fits to the IR SED. To study 
the IR emission from these two sources, we have fitted the 
SEDs using combinations of the following two models: 
(1) The inclined AGN dust torus model (Andreas Efs- 
tathiou, priv. comm.) is based on previous AGN models 



(Efstat hiou, H ough & Young 1995; Efstathiou & Rowan- 



Robinson 1995). The model contains 15 different line of sight 



orientations, with a torus opening angle of 60° and a UV op- 
tical depth of tuv = 1000. 



(2) Th e pure starburst model (Efstathiou, Rowan- Robinson 



& Siebi mmorgen 2000 ) is based on an ensemble of HII re- 
gions including an evolving population of young stars. It in- 
cludes a simple model for the evolution of HII regions, and a 
dust grai n model incorporating Polycycl ic Aromatic Hydro- 
carbons ( Biebenmorgen fc Kruegel 1992 ) . The set of models 



used have a star formation rate exponentially decaying at a 
fixed rate (r = 20Myr), but with varying initial UV optical 
depth (4 discrete values) and starburst lifetime (7 discrete 
values) . 

For PG1148+549 we have compiled ad dition al IR pho- 
tometry to that used in Rowan- Robinson (2000). We com- 
piled the SED using IRAS fluxes, together with pho tometry 
from EUV spectra ( |Tripp, Bechtold fc Green 1994| ). The / 
band flux measured in this paper is in good agreement with 
this EUVE photometry. PG1634+706, although previously 
modelled in the IR using a suite of greybodies, has not been 
previously modelled using the latest generation of radiative 
transfer codes. For P G1634+706 the F IR emission has been 
well sampled by ISO (Haas et al 1995), and we also include 
the IRAS fluxes. For both sources we also use the / band 
fluxes taken from the HST images in this paper. 

To fit the SEDs we have used the two models described 
above, allowing each model to vary freely in normalization. 
The starburst model is also allowed to vary in the discrete 
values of initial UV optical depth and starburst age. The 
torus model is also allowed to vary in inclination angle. The 
two models are then added together, and the reduced \ 2 
estimator is calculated for the total of every combination 
of starburst plus AGN model. The best fit model is then 
selected using the miniumum \ 2 - The SEDs and best model 
fits for these two sources can be found in Figure ^. 

For PG1148+549 the fit to the IR SED is excellent. 
For PG1634+706 the fit is again very good, except that it 
underpredicts the / band flux and very slightly overpredicts 
the FIR emission. The small inconsistencies in the far-IR is 
likely due to having only a small number of input models 
with discrete parameter values, which do not sample the 
full parameter space of either the starburst or AGN models. 
The underprediction of the / band flux in the SED is most 
likely also due to this, but could also be accounted for by an 
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F00235+1024, F10026+4949 and F15307+3252 are included here 
for completeness, in these cases the 'host' is in fact the source 
itself. All reliably determined profiles were ellipticals. Although 
a host for F14218+3845 was resolved the profile was undeter- 
mined. The morphology of F00235+1024 is too disturbed to fit 
a profile to. 2a errors are quoted. 1 k correction for the F814W 
filter. 2 Absolute magnitudes including k corrections. 

unquantified contribution from evolved stellar populations 
which are not included in either model. 



4.3 Companion Sources 

An L* galaxy at the highest redshift in our sample (z = 1.33) 
would have a magnitude of approximately mi = 22.2, and 
the faintest galaxies visible (to a 3<r detection) in any of our 
observations have a magnitude of mi ~ 22.5. Companion 
sources, if present, should certainly be visible for all the 
sources in the sample. 

Eight of the nine objects in the sample show at least 
one other object within the PC field of view. Two objects, 
F12358+1807 and F12509+3122 show a single large, bright 
companion whose magnitude and separation from the source 
plausibly suggest a cluster environment. F12358+1807 also 
has a number of smaller, dimmer companions. The other six 
objects all show at least two small, dim objects within 15". 
For PG1634+706 there appear to be no other sources visible 
within the PC field of view. At the redshift of this source 
(z — 1.33), this would mean that there are no companions 
with an optical luminosity comparable to or brighter than 
L* within 0.5Mpc, suggesting that this source is isolated. 



4.4 QSO Host Galaxy Properties 

Six objects in our sample are QSOs as determined from the 
optical morphologies, and host galaxies are clearly resolved 
in four of these six objects. Images of the four host galaxies 
can be found in Figure |^. For a positive host galaxy detection 
our original imposition of a minimum 95% confidence limit 
was found to be easily satisfied, all detections achieving at 
least a 97.5% confidence with two exceeding 99% confidence. 
For the two PG quasars a combination of distance and ex- 
tremely luminous central point source meant that the host 
was unresolved. The average value of the QSO host galaxy 
magnitude was found to be: 



(Mi) 



-24.9 ±0.4 



(4) 



which is comparable in value to the magnitudes of the 
brightest cluster elliptical galaxies found in the Virgo and 
Coma clusters ( Caon, Capaccioli fc D'Onofrio 1994 ; Jor 



gensen & Franx 1994). It is also approximately 1.7 magni- 
tudes brighter than the field galaxy magnitude, Mj = —23.2 
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Figure 3. IR Spectral Energy Distributions for PG1148+549 & PG1634+706. The filled triangles are IRAS fluxes, and the unfilled 
triangles with an arrow below them are IRAS upper limits. Filled diamonds are the / band fluxes from this paper. Unfilled diamonds are 
any additional data taken from the literature. For PG1148+549 the additional data is taken from Tripp et al (1994). For PG1634+706 
the additional data is taken from Haas et al (1998). 



(Pozzetti, Bruzual fe Zamorani 1996; Efstathiou, Ellis & Pe- 



terson 498q ). For comparison the three objects which do not 
possess a bright central source are, on average, about 0.5 
magnitudes brighter than the mean QSO host galaxy mag- 
nitude, with (Mj) = -25.5 ± 0.9. 

For three of the four QSOs in which a host was de- 
tected, the residual light profile was sufficiently bright and 
extended to allow a galaxy profile to be reliably fitted. All 
three hosts were found to be very luminous ellipticals, with 
large derived scalelengths. These values can be found in Ta- 
ble]! The scalelength for LBQS1220+0939, at ~ 88Kpc is 
extremely large, and is comparable to the largest observed 
ellipticals. Measured ellipticities were found to be very small, 
generally 0.1 < e < 0.2, justifying our use of a fixed elliptic- 
ity in extracting surface brightness profiles. 



4.5 Individual Sources 

4-5.1 F00235+1024 

z = 0.58. There is no X-ray detection down to a limit of 
Lx/Ltoi = 2.3 x 10 -4 , indicating either atypically weak 
X-ray emission or an obscuring column density of Nh > 



10 cm ( Wllman et al 1998). The IR power source is pre- 



dominantly of starb urst origin, but w ith a significant (37%) 
AGN contribution (IVerma et al 200 1|). 



The HST image shows an object which is clearly 
strongly interacting. There is a compact, bright region to 
the north, with traces of nebulosity extending further to the 
north. To the south the emission is less bright, more ex- 
tended and shows some small plumes. Present throughout 
this region is a series of very compact, bright, 'knots'. There 
are two small companion sources; one 6.9" to the north and 
one 23" to the east. 



4.5.2 F 10026 +49 49 

z = 1.12. The IR power source of this object is best ex- 
plained as being AGN dominated, but with an uncertain 



starburst contribution ( Rowan- Robinson 2000| ). The 'warm' 
IR colour (/25//60 = 0.665) also implies AGN dominance. 
The object is small (less than 1" in diameter), has a very low 
ellipticity (~ 0.1) and has no bright central source. There is 
some morphological disturbance in the form of faint nebu- 
losity extending towards the east. Six small, dim companion 
sources are visible within a 6" radius, two of which lie within 
3". These very close companions and asymmetric light dis- 
tribution imply recent or ongoing interactions. The surface 
brightness profile is best fitted by an elliptical profile, with 
Mi — —25.5 and a scalelength of 6.7 Kpc. 



4.5.3 PG1148+549 



z = 0.97. This radio quiet QSO is cited (Hamann et al 199 
as probably containing the broad emission line Ne VIII. This 
would require a total column of Nh > W 22 cm~ 2 , and at 
least 33% obscuration of the central regions. The presence 
of this line would mean that gas in the broad emission line 
regions would be extremely hot, and that this component 
would appear as an X-ray warm absorber if it lies along 
the line of sight to the X-ray continuum source. Long term 
broad b and optical photometry of this object (Barbieri et 
al 1988 ) shows possible evidence for a variable source, with 



suspicious small amplitude variability in the period 1968 - 
1986. 

The HST image shows a bright QSO with some faint 
nebulosity extending approximately 2.2" to the north east. 
Three companion sources can be seen 4.3" to the south 
west, 12.2" to the south and 13.2" to the south east. A host 
galaxy was detected at greater than 97.5% confidence, how- 
ever fitting a galaxy profile to the underlying light profile 
proved impossible. The measured host magnitude was also 
extremely large, at Mi — —27.3. It is therefore more likely 
that we are merely detecting a profile of light from the AGN 
that has been scattered within the host galaxy, rather than 
the host galaxy itself, hence the host is unresolved. The IR 
SED is best fit by an AGN (70%) at 65° line of sight orienta- 
tion, and a 37Myr starburst (30%) with tuv = 660. The fit 
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Figure 4. Optical spectrum for F14218+3845, taken by Andy 
Lawrence with the FOS-I spectrograph on the Isaac Newton Tele- 
scope as part of the IRAS FSS-z survey. 



shows that the torus dominates the emission from 60/im to 
the UV. The starburst dominates the emission longward of 
60/im. The orientation of the AGN model implies that the 
Broad Line regions should be visible in the optical, which is 
consistent with the optical morphology. 



4.5.4 LBQS1220+0939 

z = 0.68. This very compact QSO possesses largely symmet- 
rical, dim, extended nebulosity and appears to be morpho- 
logically undisturbed, except for some faint inhomogeneities 
to the south west. There is a very dim companion source 
2.1" to the south west and two brighter companions < 11" 
to the north east. A PSF + host galaxy fit is preferred 
to a pure PSF fit at > 99% confidence, with an elliptical 
being strongly preferred to a spiral. The host galaxy has 
Mi = —25.0, and a scalelength of ~88Kpc. The average el- 
lipticity of the source is 0.2. There is a gradient of decreasing 
ellipticity from the central regions ((e) = 0.1 at 0.25") to 
the outer regions ((e) = 0.3 at 0.80"). There are however no 
large changes in ellipticity between adjacent isophotes and 
there is no evidence for isophotal twisting, suggesing that 
the host is undisturbed. There is insufficient IR photometry 
to draw any conclusions about the IR power source. 



4.5.5 F12358+1807 

z — 0.45. This moderately bright QSO shows no signs 
of tidal interaction, the surrounding nebulosity being ex- 
tremely symmetric. There is however a bright spiral galaxy 
8.6" to the north whose magnitude (mgi4 = 20.05) sug- 
gests that it may be a companion rather than a chance 
foreground object. The spiral appears to be morphologically 
undisturbed and the spatial separation between it and the 
source is too great for direct interactions. The IR power 
source of this object cannot be constrained as there is insuf- 
ficient IR photometry available. 



The host galaxy is formally detected at > 99% confi- 
dence and an elliptical profile is strongly preferred to a spiral 
profile (x 2 ~ 1.0 vs. \ 2 ~ 4.0). Fitting elliptical isophotes 
shows that the ellipticity is of the order 0.1 - 0.2. The host 
has an absolute k corrected magnitude of -24.5, and a de- 
rived scalelength of 7.1Kpc. 



4.5.6 F12509+3122 

z = 0.78. This QSO is very symmetrical and shows no signs 
of interaction. There is however a large, bright companion 
galaxy 15.5" to the south whose magnitude (mi = 17.9) 
suggests physical association with the source rather than 
chance spatial coincidence. The companion appears to be 
morphologically undisturbed, but has a bright nucleus. The 
light distribution of the companion is best fit by an ellipti- 
cal profile, the HST image suggests it may be an SO galaxy. 
Interestingly, this galaxy, if it is a companion, has a compa- 
rable magnitude (Mi = —25.2) to the QSO host galaxies, 
and is also similar to the brightest ellipticals found in mod- 
erate to rich galaxy clusters. In addition to the bright com- 
panion there are also four other smaller, dimmer companion 
sources. 

The host galaxy is formally detected with > 97.5% con- 
fidence, has an absolute k corrected magnitude of -25.2, and 
a scalelength of 10.3Kpc. The 60/xm-25/xm colour suggests 
that an AGN interpretati on of the IR SEP is m ost likely, 
although previous studies ( Rowan- Robinson 200C ) imply an 
unquantified, but significant starburst component. The mea- 
sured ellipticity of the host was found to be small, (e) = 0.15. 



4.5.7 F14218+3845 

z = 1.21. Overall, the IR emission from this object is best 
expl ained as being pre dominantly (74%) starburst in ori- 
gin (Verma et al 2001). The starburst dominates only at 
wavelengths longward of 30/xm. Below 30/im the AGN torus 
dominates. The fit predicts that an AGN should be responsi- 
ble for most of the optical emission, which is consistent with 
the optical morphology. The optical spect rum of this source , 
taken as part of the IRAS FSS-z survey (Oliver et al 199£ ) 
on the Isaac Newton Telescope using the Faint Object Spec- 
trograph (FOS-I), is given in Figure ^. The single broad line 
is most plausibly interpreted as Mgll. Another possibility 
is H/3, but the absence of an Ha line at longer wavelengths 
makes this interpretation unlikely. There are two small dim 
companions 5.8" away to the south east and south west, and 
one slightly more extended companion 13.5" to the south 
west. The source itself is very small. A host galaxy was de- 
tected at > 99% confidence, although it proved impossible 
to distinguish between a de Vaucoleurs or spiral galaxy pro- 
file as both profiles produced an equally good fit to the data. 
The host galaxy magnitude (Mi = —24.7) is much closer in 
value to the total source magnitude than for any of the other 
sources. There is an extremely faint inhomogeneity about 
0.5" to the southeast in the form of a very small plume of 
nebulosity, this feature is however barely detected in the 
HST image (~ 3ct) and the PSF subtracted image shows no 
signs of disturbance. We thus conclude that this object is 
not involved in ongoing interactions. 

On prima facie grounds, the HST image of this source 
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shows possible gravitational lensing signatures. Further in- 
vestigation however showed no evidence for lensing. The very 
faint 'ring' surrounding the QSO is at the same radius as the 
Airy ring in the normalised TinyTim PSF generated for this 
object, and is therefore virtually certain to be a PSF arti- 
fact. The likelihood of seeing a symmetrical Einstein ring 
with perfectly uniform luminosity around the entire circum- 
ference is in any case fairly small, and would require perfect 
or near-perfect alignment between a circularly symmetric 
lens and an extended source. There is no tangential stretch- 
ing of the source, and there are no other structures or nearby 
companions that suggest any plausible lensing scenario. This 
implies that the derived IR luminosity is intrinsic to the 
source. 



4.5.8 F15307+3252 

z = 0.93. This object was first identified as an HLIRG 
by Cutri et al fll994| ). It was, after F10214+4724 and 
P09104+4109, the third HLIRG to be discovered, and bears 
many similarities to its predec essors. Optical spectropo- 
larimetry by Hines et al (1995) shows a highly polarized 
continuum, and a broad emission line region that is typi- 
cal of QSO's. Nearby companion s have been detec ted in K 
band images (3oifer et al 1994). Liu et al (1996) discuss 



Table 3. Photometry for features within F15307+3252 



the properties of this object and, using deep Keck K band 
imaging combined with 1.1/xm - 1.4/mi spectroscopy, dis- 
cuss whether this object is gravitationally lensed or a giant 
elliptical in the process of formation . Their results are incon- 
clusive. ROSAT HRI observations (Fabian et al 1996) show 



no X-ray emission down to a limit of 2 x 10^L(, o ; indicating 
either atypically weak X-ray emission or that little of the 
X-ray flux is scattered into our line of sight. 



CO line and rest frame 650/mi observations (Yun & 



Scoville 1998 ) derive upper limits to the CO molecular gas 
mass of 5 x 10 9 Mq, and a dust mass of 0.4 — 1.5 x 10 s Mq. 
Both these values are lower than those typical for the more 
gas rich IR galaxies. F urther observatio ns of this object show 
a Seyfert 2 spectrum ( Cutri ct al 1994| ). This w as later con- 
firmed by mid IR spectroscopy ( Evans et al 1998), which also 
found a molecular gas mass comparable to or less than that 
found in the most gas rich galaxies observed locally, that the 
relatively 'warm' 60/im/100/im colours imply that most of 
the IR emission emanates from a relatively small amount of 
warm circumnuclear dust (Md ust = 10 5 — 10 7 h~ 2 Mq which 
implies a molecular gas mass M gas = 10 7 ' 3 — 10 9 ' 3 K~ 2 Mq), 
and that the narrow emission lines and Li r /Lboi = 0.90 im- 
ply that the AGN is mostly obscured. The non-detection of 
CO poses a problem fo r a purely starburst inte rpretation 
( Rowan- Robinson 2000] ), and IR SED modelling ( Verma et 
al 2001) implies that an AGN is responsible for the majority 
(68%) of the IR emission. 

The HST image shows an interacting system with three 
bright components. The most prominent source is large, 
bright, and has some small, bright arclike structures imme- 
diately to the south west. The two smaller components are 
located 1.9" to the south east, with the southerly one being 
noticeably brighter than the northern one. Extremely faint, 
diffuse emission is also present surrounding all three com- 
ponents. There is one small dim companion and one small 
bright companion 7.3" and 13.2" to the south respectively. 

The extremely high resolution image of this source 
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1 refers to the large north western source, 2 refers to the small 
south eastern source and 3 refers to the small north eastern 
source. The two companion sources clearly possess different I — K 
colours, ruling out a quadruple lensing scenario. 1 No applied k 
correction. 2 Taken from Liu et al 1996. 



taken by HST, in conjunction with the Keck K band image, 
allows the best possible chance for determining whether or 
not this system is lensed. Both the Keck image (see also fig. 
3 of Liu et al 1996) and the HST I band image are given for 
comparison in Figure ^. The two possible lensing scenarios 
given by Liu et al are: (1) The large, bright elliptical com- 
ponent is the lensed nuclear structure of a z — 0.93 galaxy, 
with the brighter of the two companion sources being the 
foreground lens. (2) This source is a quadruple image grav- 
itational lens system. By comparing the HST image to the 
Keck image (kindly provided by Michael Liu) it is possible 
to examine both these possibilities. 

We first examine the possibility of a quadruple image 
gravitational lens system. Firstly, there is no evidence from 
the / band image for the lens itself, although this in par- 
ticular does not help to rule out a lensing scenario. Table ^ 
gives the magnitudes in the / and K bands for the three com- 
ponents in the system observed in both the HST and Keck 
images. The two smaller sources in the system possess differ- 
ent I — K colours. Although differential reddening through a 
foreground lens could account for colour differences between 
the four components the 'arclike' structures to the south 
west of the large elliptical component would be expected to 
be reddened in comparison to the other three components 
in the system if this were the case. Instead they are ob- 
served to be significantly bluer as they do not appear at all 
in the K band image, but are prominent in the / band im- 
age. The large, bright elliptical component does not appear 
stretched tangentially, but shows both radial and tangential 
extension. Explaining radial extension in a quadruple image 
gravitational lens is not possible unless extremely complex 
lensing scenarios are invoked, scenarios that are not plau- 
sible from the evidence in the / band image. Although a 
quadruple lens scenario cannot be completely ruled out, it 
seems overwhelmingly unlikely. 

If the bright elliptical component of the system is as- 
sumed to be the lensed galaxy, with the brighter of the two 
companion sources assumed to be the lens, then the lack of 
counter images to the south of the lens is difficult to explain. 
The smaller of the two companion sources cannot be such a 
counter image as it appears on the same side of the lens as 
the bright elliptical, such a scenario would require unrealistic 
amounts of cosmic shear. This rules out the possibility that 
the large, bright elliptical component is the lensed nuclear 
structure of a z = 0.93 galaxy, with the brighter of the two 
companion sources being the foreground lens. We conclude 
therefore that this object is not gravitationally lensed and 
that the IR emission is intrinsic to the source. 

Fitting elliptical isophotes to the largest component in 
this source shows that the ellipticity varies 0.25 < e < 0.4 as 
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Figure 5. /fecA; A" band image of F15307+3252 (left, supplied by Michael Liu) and HST I band image. Orientation and scale are the 
same as in Figure []] 



the radius 0.2 < r < 0.32 in arcseconds. The best fit profile 
of the large bright component is found to be an elliptical, 
with an associated scalelength of 12.4Kpc and an absolute 
k corrected magnitude of Mi = —26.4. 



4.5.9 PG1634 + 706 

z — 1.33. This radio quiet object is the most IR luminous op- 
tically selected QSO in the IRAS catalogues. The extraordi- 
narily high disk luminosity (~ 10 40 IFs _1 ) requires a central 
black hole of 10 10 M n , even at the Edd ington accretion limit 



(Rachen, Mannheim & Biermann 1996) . An AGN dust torus 



model is consistent with both the ISO (Haas et al 1998) and 
IRAS data, with a limit set on any starburst contribution, a 



limit whic h is confirmed by the non -detection of CO ( Evan: 
et al 1998| ). The near IR spectrum ( Evans et al 1998 ) is not 



consistent with a Seyfert 2 galaxy, and the derived molecu- 
lar gas mass is comparable to or less than that found in the 
most gas rich galaxies observed locally. They also find that 
the relatively 'warm' 60/im/lOO/^m colours imply that most 
of the IR emission emanates from a relatively small amount 
of warm circumnuclear dust (Mdust = 10 s — 1Q 7 h~ 2 Mq, im- 



10 JJ h-^M, 



e) 



plying a molecular gas mass M gas — 10 
and that the broad emission lines and Li r / Lboi = 0.35 sug- 
gests that most of the AGN is expo sed. ASCA spectroscopy 
of this object (Nandra et al 1995) agrees with these con- 
clusions to some extent. Both the FeKa line and 'hard tail' 
are not present in the spectrum, features which are almost 
universal in Seyfert 2 galaxies. The accretion disk is thus 
concluded to be either highly ionised or optically thin, both 
of which would indicate that accretion is occurring close 
to or at the Eddington accretion limit. An optically thick 
molecular torus is concluded to be either absent or subtend- 
ing a sm all angle to the X-ray source. ROSAT P SPC obser- 
vations (Rachen, Mannheim & Biermann 1996) show this 



source to have a steep soft X-ray slope, a natural explana- 
tion of which is that the soft X-ray emission is of unsaturated 
comptonisation origin, covered by a flat nonthermal X-ray 
compone nt ab ove 1.4 ± 1.1 KeV. Barvainis, Lonsdale & An- 
tonucci (1996) confirm previous radio observations of this 



source which show a high frequency excess, attributed to a 
compact synchrotron source most commonly found in radio 
loud objects. 

The HST image shows a large, bright, symmetrical 
source with no surrounding nebulosity. There are no com- 



panions within the 36" field of view. A host galaxy was de- 
tected for this object, and a spiral disk profile is preferred 
over an elliptical profile. The host was however only detected 
at ~ 80% confidence over a pure PSF fit, and the derived 
host galaxy magnitude, at Mi = —29.4, is also unphysical. 
It is therefore more likely, as in the case of PG1148+549 that 
we are merely detecting a profile of scattered light from the 
AGN and that the host galaxy is unresolved. The IR SED is 
dominated by the AGN torus component (83%), with line of 
sight inclination of 90°. The starburst model of age 6.6Myr 
and tuv = 220 is required only to explain the emission long- 
ward of 60/im. 



5 DISCUSSION 
5.1 The Sources 

5.1.1 Morphology 

On an ad hoc basis the sample can be divided into two 
groups on the basis of their morphology, namely an 'inter- 
acting' group and a 'quasar' group. The 'interacting' group 
comprises three objects; F00235+1024 and F15307+3252 
are both strongly interacting, F10026+4949 shows slightly 
weaker signs of interaction but has numerous nearby (<5") 
companions. The 'quasar' group comprises the other six 
sources, five of which are optically bright QSOs and one 
of which, F14218+3845, is a QSO where the point source is 
relatively underluminous. 

There is a strong morphological resemblance between 
the 'interacting' group of HLIRGs and previous HST imag- 
ing of ULIRGs ( jgurace et al 199$ |Borne et al 200C| ). The 
strong resemblance between three sources in this sample 
with previously observed interacting ULIRGs, and that the 
rest of the sources are QSOs, raises the question as to 
whether any of the QSOs in the sample have arisen from 
interactions or mergers. It has been suggested (Sanders et 



al 1988) that ULIRGs represent a transition phase between 



galaxy mergers and optically luminous quasars. The forma- 
tion of quasar nuclei in ULI RGs h as been investigated by 
Taniguchi, Ikeuchi & Shioya (1999). They conclude that a 
supermassive black hole (SMBH) of mass > 10 s Mq can 
form if either one of the progenitor galaxies contains a 'seed' 
SMBH of mass > 10 7 Mq undergoing efficient Bondi type 
gas accretion during the merger, or by collapse of some star 
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Figure 6. Colour - colour plot for the 9 sources in the sample. K 
corrections have been applied to the / band and IR fluxes so that 
the plot approximates rest-frame colours. The circles are those 
sources classified as QSOs on the basis of their optical spectra, 
and the squares are the interacting sources. 



cluster containing compact remnants (black holes or neu- 
tron stars) o n a time scale ~ 10 9 years. A more recent study 
of ULIRGs ( |Farrah ct al 200 1|) using HST, proposes that 
ULIRGs as a class are not simply a transition stage between 
galaxy mergers and QSOs, but instead are a much more 
diverse population whose evolution is driven solely by the 
morphologies of the merger progenitors and the local envi- 
ronment. 

Signs of interaction in the immediate vicinity of the 
QSOs in the sample would be hard to detect due to the 
bright central source, and the relatively high redshifts of 
the objects. A significant fraction of ULIRGs in which QSO 
activity is apparent have shown clear signs of interactions 
(Leech et al 1994; Clements et al 1996|) , and observ ations of 
QSOs from the PG quasar survey (Clements 200C) showed 
that QSOs with disturbed hosts were significantly more 
likely to have a higher FIR luminosity than those QSOs with 
undisturbed hosts. These studies have however involved ob- 
jects at redshifts < 0.2. Of the six QSOs in this sample, three 
show faint signs of morphological disturbance, two are ap- 
parently undisturbed and one (PG1634+706) is too distant 
to tell. Given that all the objects that do not have an optical 
QSO show clear signs of interaction, and that the other six 
are QSOs it seems plausible that the QSOs in the sample 
could represent the last stage in a merger event. HLIRGs 
would then be plausibly interpreted as a simple extrapola- 
tion of the ULIRG population to extreme luminosities. From 
the HST and IR data alone it is however not possible to rule 
out that some or all of the QSOs are not a simple extrapola- 
tion of the ULIRG population to extreme luminosities and 
that the IR emission is not triggered by interactions. 



5.1.2 An Optical/IR Comparison 

The IR power sources for those objects with previously mod- 
elled 1 - 1000/im IR SEDs can be found in Table |. For the 
five sources with sufficiently comprehensive IR photometry 



to allow detailed modelling of the IR SED, two are star- 
burst dominated and three are AGN dominated, although 
all five require both a starburst and an AGN to explain 
the IR emission. Two sources with less comprehensive IR 
photometry which cannot be modelled using detailed SED 
models show evidence for an AGN dominating the IR emis- 
sion, with a smaller, but unquantified contribution from a 
starburst. The remaining two sources have an unknown IR 
power source as they have only one IR data point. 

Table |l| gives a direct comparison between the HST mor- 
phology and the best fit model to the IR SED. Of the two 
sources that are starburst dominated, one (F00235+1024) 
is interacting and the other (F14218+3845) is a QSO. For 
the five systems where an AGN is the dominant IR power 
source, two (F10026+4949 & F15307+3252) are strongly in- 
teracting and the rest are QSOs. It can be argued that all 
the starburst systems are triggered by interactions, one is 
in ongoing interactions, whilst the second (F14218+3845) 
shows no clear signs of interaction but does possess close 
companions. For the five AGN dominated systems however, 
two are in ongoing interactions, but three are optically lu- 
minous QSOs with no signs of interaction or very close com- 
panions. There is therefore no clear correlation between the 
IR power source and the optical morphology in this sample 
of HLIRGs. The number of sources considered here is how- 
ever too small to draw truly robust statistical constraints on 
the HLIRG population. This means that these conclusions 
are by nature quite speculative, and should be regarded with 
some caution. 

For the five sources with well sampled SEDs, the 
optical morphology correlates well with that predicted 
by the dominant contributor to the SED in the opti - 
cal. For PG1148+549, F14218+3845 flVerma et al 200 1|) 



and PG1634+706 the SED fit predicts that the optical 
emission is dominated by an AGN, which is in excellent 
agreement wit h the optical imag es. For F00235+1024 and 
F15307+3252 flVerma et al 200l| ) the SED fit predicts that 
the AGN does not contribute significantly to the optical 
emission. Although the unquantified contribution from old 
stellar populations to the / band flux means that the optical 
data cannot be used to constrain the properties of the star- 
burst, the absence of AGN both in the optical images and 
in the optical region of the SED are consistent with each 
other. 

A colour-colour plot (log(S25 1 'Seo) vs. log(Seo/Ssi4 
for the nine sources in the sample can be found in Figure 
All of the optically luminous QSOs are localised in the bot- 
tom right hand corner of the colour-colour plane. The three 
interacting sources are also localised in the log(Seo/ Ssii) 
axis, but are more spread out along the log(S25/S^o) axis 
due to their varied AGN fractions. F14218+3845 is both sub- 
stantially less luminous in the optical than the other QSOs 
in the sample (probably due to higher levels of obscuration 
in the nuclear regions), and contains a starburst compo- 
nent con tributing dominan tly to the IR emission longward 
of 60^im ( Verma et al 2001), hence its position in the top left 



of the colour-colour plane close to the interacting sources. 



5.1.3 Gravitational Tensing 

The extraordinarily high bolometric IR luminosities ob- 
served in these objects suggests that magnification via grav- 
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Figure 7. HST F814W images of the four resolved HLIRG QSO host galaxies. Residual PSF features have been masked out. 



itational lensing may play a significant role. Gravitational 
lensing has been observed previously in four IR selected hy- 
perluminous sources, these include the 'Cloverleaf quasar 
and F 102 14+4724. Signs of lensing were therefore a pri- 
mary search objective for this sample. One of the sources 
was known to be a strong candidate for lensing, as mor- 
phological information for F15307+3252 from ground based 
imagi ng suggested t hat it may be a lensed system of some 
kind (|Liu et al 199^). 



It came as quite a surprise therefore, when none of the 
sources in the sample showed any evidence for gravitational 
lensing. The previous observations of HLIRGs that have de- 
tected lensing using HST indicate that this phenomenon is 
at least present to some esxtent in the HLIRG population. 
Including our sample suggests that only a small minority 
(around 15% - 20%) of HLIRGs have been mistakenly clas- 
sified as such due to the effects of lensing. This percentage 
should be treated with a degree of caution; the sample pre- 
sented here is at a lower median redshift than that of previ- 
ously confirmed lensed sources, the sample of HST imaged 
HLIRGs is also as a whole not homogenous. There is also 
the possibility that these objects may be lensed by a fore- 
ground cluster, such a lensing scenario might be expected 
to produce weaker signs of lensing than would be observed 
if any of these objects were being lensed by a single source. 
This possibility can only be fully examined by a quantita- 



tive study of the environments of these objects, although 
the complete lack of lensing signatures in any of the objects 
makes a cluster lens scenario unlikely. 



5.2 The Host Galaxies 

Previous surveys of QSOs using HST have generally been 
extrem ely successful in detecting host galaxies. Urry et al 
( |1999| ), g ve results for a sample of BL Lacertae objects, de- 



tecting host galaxies out to a redshift of z — 0.7. They find 
that the hosts are all extremely luminous ellipticals with 
{Mi) = —24.6. The derived ellipticities are generally low 
and morphologies are smooth. Imaging of a sample of 14 
quasars (6 radio loud, 5 radio quiet, 3 radi o qui et ultra- 
luminous infrared quasars) by Boyce et al (11998) showed 
that all the radio loud quasars had elliptical hosts, the radio 
quiet quasars possessed both elliptical and spiral hosts, and 
that the radio quiet ultraluminous infrared quasars lay in 
violently interacting systems. Host galaxies were on average 
0.8 mag brighter than L* . K ban d imaging of of a sample 
of 14 luminous radio quiet QSOs (Percival et al 2001) found 



that two of the host galaxies were violently interacting, and 
the rest were undisturbed ellipticals. Results from imaging 



of two luminou s radio quiet quasars (Bahcall, Kirhakos & 



Schneider 1996) show that both quasars reside in appar- 
ently normal host galaxies, one being an elliptical and one a 
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Figure 8. A simulation in which an observed, saturated Plane- 
tary Camera / band PSF is added to the large elliptical compo- 
nent of F15307+3252 shown in Figure [j], A theoretical TinyTim 
PSF has then been subtracted off, following the procedures de- 
scribed in §3.2. Residual PSF features have been masked out. The 
weak signs of morphological disturbance seen in the original image 
are still clearly visible. The less luminous of the two companion 
sources is visible uii Lhu lowei light uf the image. 



traction. In order to make the simulations realistic we took 
the large bright elliptical component of F15307+3252, which 
shows weak signs of morphological disturbance, and added 
on an observed saturated PSF taken with the F814W filter 
from the HST PSF archive. TinyTim PSFs were then gen- 
erated for the position of this observed PSF. Subtraction 
and fitting procedures were carried out as for the QSOs in 
the sample. A host galaxy was detected at > 97.5% confi- 
dence, and an elliptical profile was clearly preferred to a spi- 
ral profile. The measured host galaxy magnitude was found 
to be Mi = —24.9 ± 0.2, with an associated scalelength of 
r e — 12.8 ± 0.8 Kpc. The la errors of both these values 
lie within the original derived values for this object. Fig- 
ure ^ shows the galaxy after PSF subtraction. The signs 
of disturbance in the original image are still clearly visible. 
We therefore conclude that, had any of the detected host 
galaxies of the QSOs shown signs of interaction of a level 
comparable to that seen in the large elliptical component of 
F15307+3252, we would have detected them. 

The derived host galaxy properties are similar to the 
properties of the hosts of BL Lac objects observed by Urry 
et al. They are also similar to the host galaxy properties 
derived for radio quiet QSOs by Boyce et al and Bahcall et 
al, although we do not find any spiral hosts. Observations 
of the brightest ellipticals in the Virg o and Coma clusters 



(|Caon, Capaccioli fc D'Onofrio 1994 ; Jorgensen & Franx 



1994| ) show that they have very similar properties to our 



spiral. Disney et al (1995) present HST imaging for a sam- 
ple of four QSOs, finding that all four have elliptical hosts. 
Although they find that the morphologies of the hosts are 
all featureless they argue that the presence of multiple very 
close companions is evidence for interactions as the trigger 
for the quasar activity. 

We resolved host galaxies in seven of the nine objects, 
with the host morphology resolved in five of those cases. 
Host properties were reliably determined up to a redshift 
of 0.78, except in the case of F14218+3845 where the rel- 
atively linderhiminoiis point source allowed determination 
of the host galaxy magnitude, but not the profile. The host 
galaxy was not resolved for the two PG quasars, probably 
because the extremely bright central source in these two ob- 
jects creates a profile of scattered light which completely 
masks the host galaxy. The morphology of F00235+1024 
was too disturbed to fit a galaxy profile to. 

The five reliably resolved profiles were all best fitted 
by an elliptical profile with at least 97.5% confidence. The 
average QSO host galaxy magnitude was found to be {Mi) = 
—24.9, more than 1.7 magnitudes brighter than L* . For those 
sources that do not display a point source the galaxies were 
brighter, with (Mi) = —25.5 ±0.9. The derived scalelengths 
were also large, all but one falling in the range 6.7 < r e < 
12.5. The scalelength for LBQS1220+0939, at ~ 88Kpc is 
extremely large, and is comparable to the largest observed 
ellipticals. The possibility exists that this value is erroneous, 
possibly due to AGN light scattered within the host galaxy 
or (less likely) a poorly subtracted PSF, however the faint, 
extended emission surrounding this source is supportive of 
a large scalelength. 

We performed simulations to quantify the efficiency of 
host galaxy detection, and to determine whether signs of 
morphological disturbance would be resolved after PSF sub- 



host galaxies. Interestingly, both IRAS quasars imaged by 
Boyce et al had host galaxies classified as interacting. If 
all the QSOs in our sample arose due to interactions and 
mergers it would not be unreasonable to expect some of the 
hosts in our sample to be interacting systems, although the 
combination of a luminous point source and relatively high 
redshift may well serve to mask any such features. The low 
measured ellipticities for the QSO hosts suggests that the 
hosts are either featureless or in the very last stages of a 
merger. As ULIRGs are thought to be mergers between two 
or more gas rich spiral galaxies, and as a merger betw een two 
such galaxies is thought to form an elliptical galaxy ( Barnes 



198£), the properties of our host galaxies are consistent with 
being ULIRGs in the last stages of a merger. 

5.3 Companion Sources 

It has long been suspected that IR emi ssion in ULIRGs is 



trigg ered by interacti ons and mergers ( Banders & Mirabel 



1996), one possibility (Borne et al 200C) is that the IR emis- 
sion may arise due to multiple mergers and interactions in 
Hickson Compact Groups. Results from the analysis of the 
QSO host galaxies suggest that HLIRGs may reside in rich 
clusters, given that the hosts are very similar to the brightest 
cluster ellipticals observed locally. All but one of the sample 
show at least one companion source within 36". Two sources 
have single, large, nearby bright companions, one of which 
is a very bright elliptical, the other a spiral. The other six all 
show at least two smaller, dimmer companions within a 16" 
radius. The one exception, PG1634+706, apparently has no 
companions within 0.5Mpc, suggesting that this source does 
not reside in a cluster. 

These results, even when taken in conjunction with the 
results for the host galaxies, do not provide conclusive evi- 
dence for or against a cluster environment for all HLIRGs, 
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but do suggest that some HLIRGs may reside in a rich en- 
vironment. The environments of the HLIRGs in this sample 
will be presented in more detail in a future p aper. 



Studies of the e nvironments of ULIRGs ( Borne & Scott 



1991; Banders 1992) have shown that these objects do not 
appear to reside in regions of higher than average density. 
Most ULIRGs appear either isolated or in systems similar 
to the local group. Results from the IRAS surveys have 
shown that most ULIRGs are strong interactions between 
gas rich spirals with a relative mean velocity of < 200km 
s -1 . These are not typical conditions of rich cluster envi- 
ronments, where relative mean velocities are much higher 
and many galaxies are eit her gas poor spirals, or ellipticals 
( Sanders fc Mirabel 1991: ). A study of the cluster environ- 
ments of radio loud and ra dio qu iet AGN has been per- 
formed by Mclure & Dunlop (|200l| ). They showed that both 
classes of object resided in rich environments of approxi- 
mately Abell class 0, although there was a large scatter. 
They also found no evidence for an epoch dependent change 
in the environments, and that AGN cluster environments as 
a whole were consistent with being drawn at random from 
the general cluster population. 

The tentative hypothesis presented here that some frac- 
tion of HLIRGs may reside in a rich cluster environment, if 
confirmed, would be an important result and would raise 
problems for interpreting HLIRGs as a simple extrapola- 
tion of the IRAS galaxy population to extreme luminosities. 
Density evolution with redshift cannot explain sparse envi- 
ronments at 2 ~ 0.2 and rich environments at higher red- 
shift. In any case strong luminosity evolution, as opposed to 
density evolution, is most appropriate for expla ining faint 
radio source counts (Rowan- Robinson et al 1993). HLIRGs 



may be the cluster analogues of ULIRGs; the higher relative 
mean velocities found in rich clusters would make the merger 
driven IR luminous phase shorter, this would also provide 
a formation mechanism for cluster elliptical galaxies. This 
scenario would imply some correlation between luminosity 
and envioronment in the more luminous of the IRAS galax- 
ies. Alternatively HLIRGs may be an entirely separate class 
of object. 



6 CONCLUSIONS 

We have performed HST / band imaging for a sample of nine 
Hyperluminous Infrared Galaxies. The results are: 

(a) The morphologies in the sample include both 
strongly interacting systems, and QSOs with no clear signs 
of ongoing interaction. 

(b) No source in the sample shows any evidence for 
gravitational lensing. This includes F15307+3252, for which 
lensing had been previously suspected. Lensing is thus con- 
cluded to be entirely absent from the sample or at such a 
low level that it has a negligible impact on the observed 
luminosities. 

(c) The QSO host galaxies are all extremely luminous 
ellipticals with properties comparable to the brightest clus- 
ter elliptical galaxies observed locally. 

(d) There is no clear correlation between the IR power 
source and the optical morphologies. Of the seven sources 
with enough IR photometry to investigate the IR power 
source, two are starburst dominated and five are AGN dom- 



inated, although all seven show evidence for some level of 
AGN emission. One of the starburst dominated systems is in 
ongoing interactions, the other is a QSO with close compan- 
ions. Of the five AGN dominated systems, two are interact- 
ing and three are QSOs with no signs of interaction. There 
is however good agreement between the optical morpholo- 
gies and the dominant contributor to the optical emission 
predicted by the SED. 

(d) Two of the sources in the sample show a single bright 
companion, six show at least two small dim companions. One 
source is very probably isolated 

Previous observations of HLIRGs have seen evidence 
for amplification of intrinsic luminosities via gravitational 
lensing. None of our sources however showed any evidence 
for this phenomenon, including one for which lensing had 
been previously strongly suspected. We conclude that only 
a small minority (perhaps 15% - 20%) of HLIRGs have been 
misclassified due to gravitational lensing. 

These results are consistent with the majority of the 
HLIRG population constituting the extremely bright end of 
the ULIRG population, i.e. where interaction induced star 
formation and quasar activity lead to the very high IR lu- 
minosities. The evidence for this can be summarised as; (1) 
Three of the sources lie in interacting systems, which bear 
a strong morphological resemblance to ULIRGs (2) For the 
five sources with sufficient data to accurately model the IR 
SED, the IR emission is best explained as arising from a mix- 
ture of starburst and AGN activity, which is similar to that 
found amongst ULIRGs. (3) The QSO host galaxy proper- 
ties correspond closely to locally observed bright cluster el- 
lipticals, a cluster environment would provide a convenient 
source for interactions and mergers. (4) Most of the sources 
show at least two nearby companions. (5) One of the bright 
companion sources is a galaxy with similar properties to 
bright cluster ellipticals. 

There are however some important points to note about 
these conclusions. None of the QSOs or QSO hosts show any 
clear signs of interaction, signs which have been observed 
in ULIRGs in which QSO activity is apparent. Such signs 
may however be hidden due to high redshift, luminous cen- 
tral source, or because hyperluminous activity takes place 
in mid to late stage mergers where the signs of merging are 
diminished. The evidence that the QSOs reside in luminous 
elliptical hosts is only suggestive that some may also reside 
in rich cluster environments, and is not conclusive. For ex- 
ample one source (PG1634+706) is apparently isolated. Fur- 
ther analysis is required to make more concrete statements 
about the local environments of HLIRGs. If the tentative 
suggestion put forward here that some HLIRGs may reside 
in a rich cluster environment is found to be correct, then 
this would pose problems for HLIRGs being a simple ex- 
trapolation of ULIRGs to extreme luminosities as ULIRGs 
are known not to reside in rich clusters. It is thus possible 
that some HLIRGs are the cluster analogues of ULIRGs, or 
that the IR emission in some HLIRGs may not be triggered 
by interactions. 

It must be emphasized that these conclusions are based 
on a relatively small sample of objects, essentially those 
HLIRGs that have been observed at sufficient resolution to 
reliably determine morphologies, host galaxy properties and 
lensing characteristics. It is not possible to draw statistically 
significant conclusions from a sample of this size. Observa- 
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tions of a larger number of HLIRGs to the highest possible 
depth and resolution are the only way to confirm or refute 
the conclusions drawn above. 
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